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ABSTRACT 


This  thesis  presents  a  comparison  of  three  Magnetic 
Anomaly  Detection  (MAD)  models:  a  cro s s -cor re  1  a t ion 
detection  model,  a  square  law  detection  model,  and  a  model 
referred  to  as  the  OPTEVFOR  detection  model.  FORTRAN  and 
BASIC  programs  for  the  three  detection  models  are  included 
in  this  thesis.  The  programs  yield  detection  probabilities 
for  straight  line  encounters.  Magnetic  signal  values  for 
the  straight  line  encounters  are  an  additional  output. 
Plots  of  lateral  range  curves  and  magnetic  signal  values  are 
presented.  A  discussion  of  the  required  parameters  is 
included  in  the  thesis  to  facilitate  the  use  of  the 
programs.  The  parameters  that  were  considered  in  the 
comparison  of  the  three  detection  models  are:  magnetic 
noise,  aircraft  and  submarine  headings,  submarine 
displacement,  and  the  vertical  separation  between  submarine 
and    aircraft. 
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I.   INTRODUCTION 

This  thesis  presents  a  comparison  of  three  Magnetic 
Anomaly  Detection  (MAD)  models.  The  comparison  is  in  terms 
of  probabilities  of  detection  that  were  computed  using  the 
models.  Two  of  the  models,  the  cross-correlation  model  and 
the  square  law  model,  have  been  used  to  model  sonar 
detection  [Ref.  1:  pp.  343-357].  The  third  model,  referred 
to  as  the  OPTEVFOR  model,  is  a  slant  range  threshold 
detection  model.  The  results  of  the  comparisons  are 
presented  in  graphical  and  tabular  form.  In  addition,  plots 
of  magnetic  signals  for  selected  lateral  ranges  and  noise 
levels  are  shown.  The  effects  of  noise,  aircraft  and 
submarine  headings,  submarine  displacement,  and  vertical 
separation   are    also    indicated. 

The  models  were  implemented  using  the  FORTRAN  and  BASIC 
programs  that  are  listed  in  Appendix  A.  For  those 
interested  in  using  the  programs  for  other  investigations, 
an  input  parameter  discussion  is  provided  in  Chapter  3.  To 
use  the  FORTRAN  program,  the  user  specifies  the  input 
parameters  in  an  input  file.  After  execution  of  the 
program,      an     output     file     is     generated     that     contains 


The    programs    are    based   on  an   unpublished     BASIC   program 
by    R.N.     Forrest       for    an    H.P.-    85    microcomputer. 


probabilities  of  detection  for  each  of  the  three  models.  In 
addition,  magnetic  signal  values  and  magnetic  signal  values 
plus  random  magnetic  noise  values  for  one  of  the  encounters 
generated  by  the  program  are  included  in  the  output  file. 
An  IBM  GRAFSTAT  graphical  package  was  used  to  produce  the 
graphics  in  this  thesis. 

To  use  the  BASIC  program,  the  user  must  interactively 
enter  the  input  parameter  values  for  each  encounter.  After 
execution  of  the  program,  an  optional  hardcopy  printout 
supplies  the  input  parameter  values  and  a  table  of  detection 
probabilities  for  each  of  the  three  models  (see  Appendix  A). 
Following  this,  lateral  range  curves  are  displayed  to  the 
user  for  immediate  observation.  A  typical  program  run 
producing  21  detection  probabilities  for  each  model  requires 
approximately  10  minutes  of  computing  time  on  an  Atari  800 
microcomputer. 
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II.       MODEL    DESCRIPTIONS 

A.       THE    CROSS-CORRELATION    AND    SQUARE    LAW    DETECTION    MODELS 

The  cross -cor re  la tion  and  square  law  detection  models 
are  described  in  detail  by  Forrest  [Ref.  2:  pp.  33-35].  The 
models  are  based  on  the  following  assumptions:  the  noise  is 
gaussian,  and  the  signal  sample  points  are  such  that 
adjacent    magnetic    noise    samples    are    independent. 

The  magnetic  signal  values,  as  measured  by  an 
aircraft's  magnetometer  for  the  cross -cor re  la ti on  and 
square  law  detection  models,  are  the  submarine  magnetic 
field  values  in  the  direction  of  the  earth's  magnetic  field 
at  the  positions  of  the  magnetometer.  The  submarine's  field 
is  assumed  to  be  a  dipole  field,  and  the  aircraft  and  the 
submarine  are  assumed  to  keep  constant  speeds  and  headings 
during    an    encounter. 

For  the  cr os s- corr elat ion  model,  a  complete  prior 
knowledge  of  the  magnetic  signal  is  required. 
Operationally,  this  suggests  that  a  signal  file,  which 
contains  a  replica  of  the  signal  for  each  possible  encounter 
situation,  would  be  required.  The  model  describes  a 
perfect  detection  system  with  respect  to  the  noise  model 
that  is  used.  For  the  square  law  model,  a  signal  replica  is 
not   required.       This    model    might   be   considered    to   describe 


1  1 


the  limiting  detection  capability  for  an  automatic  system 
that  does  not  use  information  about  the  shape  of  the 
magnetic    signal. 


B.       OPTEVFOR    MAD    DETECTION    MODEL 

The  OPTEVFOR  model  is  described  by  Forrest  [Ref.  3:  pp. 
7-8].  In  characterizing  the  submarine  magnetic  signal  as  a 
simple  dipole  signal,  the  U.S.  National  Defense  Research 
Committee,  [Ref.  4:  p.  20],  reports  that  the  magnetic  signal 
of  the  submarine  "varies  as  the  inverse  cube  of  the  distance 
from  the  source".  In  an  OPTEVFOR  report  [Ref.  5:  p.  1,  encl. 
1],  the  results  of  a  regression  analysis  on  empirical  peak 
to  peak  signal  output  against  slant  range  between  submarines 
and  aircraft  are  reported.  These  results  also  suggested 
this  inverse  cube  relationship  for  the  magnetic  signal. 
This  relationship  is  the  basis  for  the  OPTEVFOR  detection 
model. 

The  model  has  a  deterministic  mode  and  a  stochastic 
mode,  each  of  which  involves  the  following  parameters:  the 
submarine  magnetic  moment  (M),  an  Operator  Recognition 
Factor  (ORF),  the  average  peak  to  peak  magnetic  noise  (N)  in 
the  operating  area,    and  a   slant  range   (R).  The  relationship 


between    these    quantities    is    given    by: 

1/3 
R= 


c    M 
(ORF)    N 


(eqn    2.1) 
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The  value  of  the  constant  c  is  0.10  for  M  in  oersted 
centimeters    ,    R    in   meters,    and   N    in   gamma. 

In  the  deterministic  mode,  detection  occurs  if  and  only 
if  the  aircraft's  slant  range  from  the  submarine  at  CPA  is 
less  than  or  equal  to  R.  This  mode  yields  a  rectangular 
("cookie  cutter")  lateral  range  curve  with  the  probability 
of  detection  equal  to  1  for  an  encounter  where  the  slant 
range  at  CPA  is  less  than  or  equal  to  R,  and  0  when  it  is 
greater    than    R. 

The    stochastic    mode   allows    a    more    uncertain   approach    to 

detection    by    allowing    a    gradual     rise     in    probability    of 

detection   as    the    slant    range    at    CPA   decreases.       In   this    mode 

one    sets     the    probability    of    detection    at    R    equal     to     50 

percent,     and    the    lateral    range    curve    is    given    by    p^    =    (J)(x); 

where     it     is     understood     that    v   is     the     standard     normal 

cumulative   distribution    function   and    x    is    determined    by    the 

following    equation: 

R    -    CPA 

x    = 

(AL)    R  .  (eqn    2.2) 

In    this    equation,    CPA    is    the    magnitude    of    the    slant    range 

distance    at    CPA,    and    R    is    the    calculated    range    from    Equation 

2.1.       The    product     (AL)R    represents    a    standard    deviation. 

The    value   of    AL    can    be    considered    to    be    determined   by      "the 

combined    uncertainty    and   variability    in    the    values    of    M,     N, 

and   ORF"    [Ref.    3:    p.    8].    Two   values    of    AL    (.20    and    .01)    are 

shown  in  Figure   2.1.      If   empirical    data   was  available,    the 
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value  of  AL  could  be  chosen  to  provide  a  best  fit  to  the 
observed  results.  Note,  as  AL  approaches  0,  the  stochastic 
mode  approximates  the  deterministic  mode. 
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Figure  2.1.   Lateral  Range  Curves  for  Different  Values 
of  AL. 
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III.       INPUT    PARAMETERS 

The  input  parameters  for  the  FORTRAN  program  are  all 
contained  in  one  input  file.  This  allows  parameter  values 
to  be  easily  changed  without  recompiling  the  main  program  or 
subroutines.  Also,  with  a  few  changes,  this  program  could 
be  altered  to  operate  in  conjunction  with  a  larger  program 
to  yield  a  probability  of  detection  on  an  individual  MAD 
run. 

The  input  parameters  are  divided  into  four  areas  for 
discussion.  They  are:  (1)  sample  interval,  (2)  earth 
magnetic  field,  (3)  submarine  moments,  and  (4)  other 
inputs . 

A.   SAMPLE  INTERVAL 

The  choice  of  a  sample  interval  is  discussed  by  Forrest 
[Ref.  2:  pp.  27-30].  In  the  program,  the  total  observation 
time  in  seconds  over  which  the  samples  are  taken  is  entered 
in  T7.  This  time  should  be  long  enough  to  encompass  a 
"complete  signal"  at  the  maximum  expected  detection  slant 
range. 

As  the  slant  range  from  the  submarine  to  the 
magnetometer  increases,  the  distance  over  which  a 
significant  magnetic  signal  is  present  at  the  magnetometer 
also    increases.      Figure    3.1    graphically    shows    the    difference 
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in  the  amount  of  time  that  a  signal  is  present  for  slant 
ranges  of  200  meters  and  805  meters.  In  this  thesis,  the 
total  time  for  a  straight  line  encounter  is  assumed  to  be  20 
seconds.  As  can  be  be  seen  from  Figure  3.1,  a  20  second 
interval  adequately  covers  the  significant  portion  of  the 
magnetic  signal  for  an  805  mmeter  slant  range  at  CPA. 


200  METERS  CPA  SLANT  RANGE 


805  METERS  CPA  SLANT  RANGE 
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Figure    3.1.       Magnetic    Signals    for    Slant  Ranges    of   200    Meters 
and    805    Meters. 

The  time  between  samples  is  set  equal  to  the  reciprocal 
of  twice  the  upper  bandpass  filter  frequency  of  the  MAD 
sensor.  A  value  of  0.9  Hz  was  suggested  for  use  by  Texas 
Instruments  [Ref.  6:  p.  112]  as  an  upper  bandpass  filter 
limit  in  a  discussion  on  the  effects  of  noise  on  a  MAD 
system.  This  value  yields  a  time  interval  between  samples 
of    0.55   seconds. 

The  sample  interval  length  and  the  false  alarm  rate  (the 
expected    number    of    false    alarms    per    hour)    determine    the 
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false  alarm  probability.  The  false  alarm  rate  (F2)  is 
assigned  a  value  of  3  based  on  a  report  by  OPTEVFOR  [Ref.  5: 
p.     2.1]. 

B.       EARTH    MAGNETIC    FIELD 

Input  values  for  the  earth  magnetic  field  intensity  and 
inclination,  or  dip  angle,  may  be  taken  from  two  Defense 
Mapping  Agency  Hydrographic  Center  charts,  [Refs.  7  and  8 
respectively],  or  approximated  by  using  a  program.  If  chart 
values  are  entered,  the  earth  field  intensity  must  be  in 
units  of  gamma  and  the  inclination  in  decimal  degrees.  The 
program  used  to  determine  the  intensity  of  the  earth  field 
and  inclination  is  based  on  a  simple  dipole  field  model  that 
is    described    by    Forrest    [Ref.     9:     pp.     39-43]. 

Table  III-1  displays  the  program  output  values  of 
inclination  in  decimal  degrees  and  earth  magnetic  field  in 
gamma  for  selected  geographic  locations.  In  addition, 
corresponding  values  obtained  from  the  Defense  Mapping 
Agency  Hydrographic  Center  Charts  Number  30  and  Number  39 
are  also  displayed.  The  last  three  columns  are  the  average 
slant  range  in  meters  at  which  a  50  percent  probability  of 
detection  is  obtained  for  the  three  program  detection 
models.  The  program  input  parameters  for  these  slant  ranges 
were  the  same  as  the  base  case,  except  for  the  following 
differences:  a  sample  interval  time  of  40  seconds,  aircraft 
and     submarine     headings     of     0     degrees,      and     a     submarine 
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displacement  of  7,000  tons.  The  correlation  between  the 
slant  ranges,  comparing  the  chart  values  and  model  values, 
was  found  to  be  95  to  96  percent  for  the  three  models.  This 
suggests  that,  even  though  differences  exist  between  the 
chart  values  and  model  values,  there  is  a  high  degree  of 
correlation    in    the    final    output. 

A  limitation  to  the  simple  dipole  field  model  is  that 
it  does  not  give  an  angle  of  declination  (variation)  with 
sufficient  accuracy.  As  a  result,  all  headings  entered 
into  this  program  must  be  in  magnetic  degrees.  The  Phoenix 
Corporation  [Ref.  10:  pp.  24-25]  reports  on  geomagnetic 
field  models  that  can  represent  the  earth  field  "with 
overall  accuracies  better  than  approximately  150-200  gammas 
in  magnitude  and  .2  in  direction  of  the  field.1  This 
degree  of  accuracy  is  not  needed  for  this  program,  but  a 
simplified  version  of  one  of  these  models  that  provided 
satisfactory  angles  of  declination  would  be  beneficial  if 
the  program  were  to  be  incorporated  into  a  larger  model  that 
utilized    true    headings    as    inputs. 

C.       SUBMARINE    MAGNETIC    DIPOLE    MOMENT 

If  a  submarine's  magnetic  dipole  moment  is  known  for  the 
geographical    location    and    the    submarine's    magnetic    heading, 


Private  communication  from  R.N.  Forrest,  who 
investigated  the  use  of  the  simple  dipole  model  for  this 
purpose . 
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the  following  values  may  be  entered  in  the  program:  (1)  P, 
its  magnitude  in  oersted  centimeters  cubed,  (2)  A,  its 
direction  in  decimal  degrees  relative  to  magnetic  north,  and 
(3)  B,  its  depression  angle  from  the  horizontal  in  decimal 
degrees.  If  it  is  not  known,  these  values  must  be 
calculated  for  a  specific  location  and  magnetic  heading.  A 
program  is  included  in  the  main  program  that  can  be  used  to 
calculate  these  values.  The  program  is  based  on  a  model 
described  by  Forrest  [Ref.  9:  pp.  35-38].  The  input  to  the 
program  is  submarine  displacement  in  tons.  The  program  also 
contains  coefficients  which  relate  displacement  to  magnetic 
moment.  The  values  used  in  the  program  are  based  on  values 
cited  by  Texas  Instruments  [Ref.  6:  p.  4]. 

The  past  history  of  the  submarine  is  represented  by  the 
permanent  longitudinal,  transverse,  and  vertical  moments  of 
the  submarine  (M4,  M5  and  M6  in  oersted  centimeters  cubed). 
For  the  examples  in  this  thesis,  it  was  assumed  that 
effective  deperming  had  been  performed  and  program  default 
values  of  zero  were  used. 

D.   OTHER  PARAMETERS 

1  .   Headings  and  Speeds 

Since  the  simple  dipole  earth  field  model  used  by 
the  program  does  not  produce  accurate  angles  of  declination, 
magnetic  headings  are  required.   In  addition,  the  headings 
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must     be     in     decimal    degrees.        The     input     parameters     for 
submarine    speed   and    aircraft    speed   are    entered    in    knots. 

2 .  Noise 

The    magnetic     noise    is     assumed    to    be    such     that 

adjacent     magnetic     noise     samples     are     independent.      This 

assumption    is    based    in    part    on    the     filtering    that     is 

performed  on    the   magnetic    signal   by    the    processing    system    in 

a   MAD   detection    sensor.      The    standard   deviation   of    the    noise 

in    units    of    gamma    is    the    value   entered    into   S1  .      This    value 

can    be   approximated    from    operational    data    by    taking    from 

one-fourth    to    one-sixth    of     the    measured    peak    to    peak 

I 
magnetic    noise.       [Ref.     2:    pp.    28-29] 

The  OPTEVFOR   detection   model    incorporates    a   value    of 

average    peak    to    peak    magnetic    noise    (N)    in    the    inverse    cube 

law     calculation.        In    the    program,     the    value     of    N     is 

determined    by   multiplying    the    S1    entry   by    four. 

3 .  Operator    Recognition    Factor    (QRF) 

The  ORF  is  the  value  of  the  ratio  of  magnetic  signal 
to  magnetic  noise  for  which  the  average  operator  would 
detect  a  signal  50%  of  the  time  in  the  presence  of 
background  noise  for  a  false  alarm  rate  of  3  per  hour.  An 
ORF  value  of  3  was  suggested  for  use  by  OPTEVFOR  [Ref.  5:  p. 
4.12]. 

4.  Distance    Parameters 

Two  parameters,  R8  and  N7,  are  used  to  define  the 
points    plotted    on    the     lateral    range    curves.        R8     is     the 
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maximum  positive  value  of  the  lateral  range  in  meters  for 
which  a  lateral  range  curve  value  is  to  be  computed.  N7 
represents  the  number  of  lateral  range  curve  values  that  are 
to  be  computed  from  the  maximum  lateral  range  to  zero 
lateral  range. 

The  vertical  separation  (Z)  is  the  sum  of  the  submarine 
depth  and  aircraft  altitude  in  meters. 
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IV.       RESULTS 

Program  outputs  of  the  three  models  for  a  set  of  base 
case  conditions  are  presented  in  this  section.  Outputs  for 
variations  from  the  base  case  are  also  presented.  The 
lateral  range  of  an  encounter  (the  horizontal  separation 
between  the  submarine  and  magnetometer  when  the  magnetometer 
is  at  CPA)  for  a  50%  probability  of  detection  is  used  as  a 
measure  of  comparison.  Signal  and  signal  plus  "noise" 
traces  for  several  cases  are  presented.  The  traces  are 
based  on  the  signal  and  noise  models  that  are  part  of  the 
cross-correlation  and  square  law  models.  These  idealized 
signal  traces  appear  to  have  the  characteristics  of  actual 
signal  traces.  This  suggests  that  the  signal  and  noise 
models    might   be    used    for    training    purposes. 

A.       BASE    CASE 

The  base  case  conditions  are  listed  in  Table  IV- 1.  The 
table  is  ordered  in  the  same  manner  that  the  values  are  read 
into  the  program.  An  annotation  of  each  entry  is  included 
for    clarity. 

Figure  4-1  presents  the  lateral  range  curves  for  the 
base  case.  Points  on  the  lateral  range  curves  are  indicated 
by  the  first  letter  of  the  name  of  the  model  from  which  they 
were  derived.      The    slight  asymmetry    of    the    cross-correlation 
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detection  model  and  square  law  detection  model  curves  is 

reflective  of  the  shape  of  the  signals  that  are  'received' 

in  these  models. 

Table  IV- 1  .   Input  Parameters  for  the  Base  Case 

1.8  twice  the  upper  bandpass  limit  in  seconds 

20.0  sampling  time  interval  in  seconds 

3.0  false  alarms  per  hour 

0  Enter  inclination  (1  =  yes,  0  =  no)? 

30.0  area  of  operation  latitude  in  decimal 

degrees 
60.0  area  of  operation  longitude  in  decimal 

degrees 
45.0  submarine  magnetic  heading  in  decimal 

degrees 
10.0  submarine  speed  in  knots 

315.0  aircraft  magnetic  heading  in  decimal 

degrees 
220.0  aircraft  speed  in  knots 

0  Enter  submarine  moment  (1  =  yes,  0  =  no)? 

0  Enter  earth  field  (1  =  yes,  0  =  no)? 

0  Enter  submarine  perm  moments  (1=  yes, 

0  =  no)? 
4000.0  submarine  displacement 

200.0  vertical  separation  in  meters 

0.1  noise  (standard  deviation)  in  gamma 

1500.0  maximum  lateral  range  in  meters 

50.0  divisions  of  lateral  range 

3.0  ORF  (Operator  Recognition  Factor) 

0.2  variability  factor  for  OPTEVFOR  model 

0025  lateral  range  iteration  number  for  the 

magnetic  signal  and  signal  plus  noise  in 

the  output  file 

Table  IV-2  lists  lateral  detection  ranges  and 
corresponding  slant  detection  ranges  at  CPA  for  a 
probability  of  detection  equal  to  50  percent  for  the  cross- 
correlation  and  square  law  detection  models.  An  equivalent 
ORF    value     for    each     model     is     also     listed.  Due     to     the 

asymmetry     of     the     lateral     range     curves     for     the     cross- 
correlation   and    square    law    models,    the    average  of    the   two    50 
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Figure  4-1.  Lateral  Range  Curves  of  the  Cross-Correlation 
(C),  Square  Law  (S),  and  OPTEVFOR  (0)  Models  for  the  Base 
Case. 


Table    IV-2.    The    Lateral    Detection    Ranges,    Slant    Detection 
Ranges,    and   ORF's    of    the   Three    Models    for    the    Base    Case. 

Lateral  Detection   Slant  Detection   ORF 
Range  (meters)      Range  (meters) 


Cross -Correlation 
Square  Law 
OPTEVFOR 


885 
685 
318 


907 
714 
376 


.21 

.44 

3 


percent  detection  ranges  was  used  as  the  lateral  detection 
range.  The  equivalent  ORF  values  for  the  cross-correlation 
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and  square  law  detection  models  were  calculated  using  the 

slant  detection  range  values   with  the  following  equation, 

which  was  obtained  from  Equation  2.1: 

c  M 

ORF  = 

R3  N   .  eqn  4.1 

For  the  base  case,  the  magnitude  of  the  submarine  field  (M) 
at  the  submarine  is  6.35  x  108  orested  cm  ,  the  noise  (N)  is 
.4  gamma,  and  the  value  of  the  constant  (c)  is  .1  .  This 
suggests  that,  in  order  to  detect  a  magnetic  signal  50 
percent  of  the  time  with  a  false  alarm  rate  of  3  per  hour, 
the  magnetic  signal  to  magnetic  noise  ratio  should  be  .21 
for  an  ideal  cross-correlation  detector  and  .44  for  an  ideal 
square    law    detector. 

Using  the  ORF  values,  the  cross -correlati on  and  square 
law  detection  models  can  be  used  to  describe  the  performance 
of  an  operator.  To  do  this,  a  modified  value  of  the 
standard  deviation  ( C )  of  the  input  noise  can  be  used.  The 
modified  value  is  equal  to  (ORF  )  (o*)  /  .2  1  for  the  cross- 
correlation  detection  model  and  (ORF)  (o*)/. 44  for  the  square 
law  detection  model.  With  these  modifications,  the  two 
models  can  be  used  to  describe  the  detection  capability  of 
an  operator  with  a  specified  ORF.  An  example  of  a  lateral 
range  curve  with  the  modified  noise  standard  deviation  for 
an  ORF  of  3  is  presented  in  Figure  4.2  for  each  model. 
These  curves  are  comparable  to  the  lateral  range  curve  for 
the   OPTEVFOR    model    that    is    shown    in    Figure    4.1. 
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CROSS-CORRELATION:  NOISE  -   1.43  GAMMA 


SQUARE  LAW;  NOISE  -  0.68  GAMMA 


Figure  4.2.  Cross-Correlation  and  Square  Law  Lateral  Range 
Curves  to  Describe  the  Performance  of  an  Operator  with  an 
ORF   of    3. 


The  automatic  MAD  system  manufactured  by  Canada's  CAE 
Electronics  Ltd.  is  expected  to  produce  a  50  percent  in- 
crease in  detection  slant  range  [Ref.  11  ].  Using  the 
detection  slant  range  for  the  OPTEVFOR  model  of  376  meters, 
a  50  percent  improvement  would  yield  a  detection  slant  range 
of  564  meters.  The  ORF  for  a  detection  system  with  this 
capability  would  be  .88.  The  cross-correlation  and  the 
square  law  detection  models  could  be  used  to  yield  lateral 
range  curves  for  a  system  with  an  ORF  of  .88  by  using  a 
noise  standard  deviation  equal  to  .88  (0")/.21  and  .88  ( 0- )  / 
.44  respectively.  Figure  4.3  shows  the  lateral  range  curves 
of  the  two  detection  models  with  a  50  percent  improvement  in 
slant  range  detection.  Note,  with  the  modified  noise 
standard  deviations,  the  models  are  essentially  equivalent 
for    the    cases    considered. 
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CROSS-CORRELATION;  NOISE  -  .2  GAMMA 


SQUARE  LAW;  NOISE  •   .2  GAMMA 
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Figure  4.3.   The  Cross -Correlation  and  Square  Law  Models  to 
Describe  LRC ' s  for  the  CAE  Automatic  Detection  System. 


SIGNAL  AT  CPA  OF  0  METERS 


SIGNAL  +  NOISE  AT  CPA  OF  0  METERS 
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SECONDS 


SECONDS 


Figure  4.4.  Magnetic  Signal  and  Magnetic  Signal  Plus 
Magnetic  Noise  at  a  Lateral  Range  at  CPA  of  0  Meters  for  the 
Base    Case. 


Figures  4.4  and  4.5  present  the  magnetic  signal  and  a 
representation  of  magnetic  signal  plus  magnetic  noise  that 
would  be  received  under  the  base  conditions  by  a 
magnetometer    with   a    lateral    range   of    0    meters    and    of    780 
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meters.  The  signal  plus  noise  trace  was  generated  from 
signal  plus  noise  values  obtained  by  adding  a  signal  value 
to  a  gaussian  noise  value.  The  gaussian  noise  value  was 
generated  by  multiplying  the  standard  deviation  of  the  input 
noise  by  a  pseudo  normal  random  number  from  a  population 
with  mean  0  and  variance  1.  The  pseudo  normal  random  numbers 
were  generated  using  LLRANDOMII,  a  resident  program  at  the 
Naval    Postgraduate    School    computer     [Ref.     12:     p.     2.2]. 


SIGNAL  AT  CPA  OF  780  METERS 


SICNA1  +  NOISE  AT  CPA  OF  780  METERS 


SECONDS 


SECONDS 


Figure  4.5.  Magnetic  Signal  and  Magnetic  Signal  Plus 
Magnetic  Noise  at  a  Lateral  Range  at  CPA  of  780  Meters  for 
the    Base    Case. 


The  magnitude  of  the  magnetic  signal  shown  in  Figure  4.4 
is  very  large  in  comparison  to  the  background  noise.  The 
peak  to  peak  signal  to  noise  ratio  is  approximately  14  to  1 . 
An  operator  would  have  little  difficulty  identifying  the 
signal    in    this    signal    plus    noise    trace. 

Conversely,  the  magnetic  signal  shown  in  Figure  4.5  is 
small    compared    to    the    background    noise.      The  peak   to  peak 
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signal  to  noise  ratio  is  .35.  The  probabilities  of 
detection  for  the  lateral  range  of  780  meters  are:  .95  for 
the  cros s- correlat ion  detection  model,  .28  for  the  square 
law  detection  model,  and  0  for  the  OPTEVFOR  detection  model. 
It  seems  apparent  that  an  operator  would  have  a  difficult, 
if  not  impossible,  time  in  detecting  this  signal  at  a 
reasonable    false    alarm    rate. 

B.   DIFFERENT  NOISE  INPUTS 

The  first  variation  on  the  base  case  shows  the  effect  of 
different  noise  inputs.  The  standard  deviation  (  O' )  of  the 
peak  to  peak  noise  is  the  input  parameter  that  is  varied. 
Table  IV-3  lists  the  different  &*  values  and  the 
corresponding  lateral  detection  ranges. 

Table  IV-3.   The  Effect  of  Noise  on  Detection  Range. 

Standard  Deviation      Lateral  Detection  Range  in  Meters 
of  Noise  in 

Gamma  Cross- 

Correlation    Square  Law      OPTEVFOR 

.005  2250(2259)*  1792(1803)*  1000  (1020)* 

.01  1832  (1843)  1446  (1460)  782   (807) 

.05  1110  (1128)      868   (890)  427   (472) 

.1               885   (907)      685   (714)  318   (375) 

.5               512   (550)      382   (431)  90   (219) 

it 

The     numbers     in     parentheses     are     the     slant     range 
distances    in   meters.      The   vertical    separation    is    200    meters. 

Figure  4.6  displays  lateral  range  curves  for  the  three 
models  when  the  standard  deviation  of  the  noise  is  .01 
gamma.        These    three    curves     show    an     increase     in     lateral 
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detection  range  over  the  base  case.  Note  that  the  asymmetry 
of  the  cross-correlation  and  square  law  detection  models  is 
more  apparent  in  Figure  4.6  than  it  was  in  Figure  4.1. 

Figure  4.7  displays  the  magnetic  signal  (which  is  the 
same  as  the  signal  in  Figure  4.5)  and  the  magnetic  signal 
plus  magnetic  noise  at  a  horizontal  distance  of  780  meters 
when  the  magnetometer  is  at  CPA.  The  signal  to  noise  ratio 
is  3.5.  The  figure  suggests  that  a  MAD  operator,  in  this 
case,  should  have  the  ability  to  detect  a  signal  at  780 
meters  lateral  range  with  a  satisfactory  false  alarm  rate. 


-2000 


2000 


Figure  4.6.   Lateral  Range  Curves  for  the  Three  Models  with 
the  Standard  Deviation  of  the  Noise  Set  to  .01  Gamma. 
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Figure  4.7.  Magnetic  Signal  and  Magnetic  Signal  Plus 
Magnetic  Noise  with  the  Standard  Deviation  of  Noise  =  .01 
Gamma    at    780    Meters  Lateral    Range. 


C.       DIFFERENT    HEADINGS 

The  headings  of  a  submarine  and  an  aircraft  in  an 
encounter  have  an  effect  on  detection  ranges.  The  effect  of 
different  headings  was  investigated  using  the  square  law 
detection  model,  and  the  results  in  terms  of  lateral 
detection  ranges  are  presented  in  Table  IV-4.  This  table 
suggests  that  a  submarine  should  choose  a  magnetic  heading 
of  either  East  or  West,  and,  for  an  encounter,  an  aircraft 
should   also   choose   a    magnetic   heading   of    East    or   West. 

Figure  4.8  shows  lateral  range  curves  for  a  submarine 
heading  North  and  an  aircraft  heading  East.  In  this  case, 
both  the  cross -correla tion  and  square  law  detection  model 
lateral  range  curves  display  noticeable  asymmetry.  The 
OPTEVFOR   detection   model    lateral   range  curve   is   symmetric 
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Table    IV-4.     Square     Law     Lateral     Detection    Ranges 
Different  Submarine    and  Aircraft    Magnetic    Headings 


for 


Aircraft 
Headings 
(magnetic ) 


0 
45 
90 
135 
180 
225 
270 
315 


0 


700 
724 
730 
712 
685 
712 
730 
724 


45 


650 
637 
654 
685 
640 
624 
646 
685 


Submarine  Headings  (magnetic) 


90 


135 


180 


225 


498 
505 
527 
505 
498 
501 
51  9 
501 


640 
670 
654 
636 
650 
672 
646 
624 


686 
713 
730 
724 
700 
724 
730 
713 


640 
624 
646 
682 
650 
636 
654 
682 


270 


498 
501 
524 
502 
498 
506 
533 
505 


31  5 


650 
676 
646 
624 
640 
672 
654 
636 


Figure  4.8.  Lateral  Range  Curves  for  the  Submarine  Heading 
North  and  the  Aircraft  Heading  East. 
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but,  like  the  curves  for  the  other  models,  it  shows  an 
increase  in  detection  ranges  over  those  for  the  base  case 
(where  the  submarine  is  heading  NE  and  the  aircraft  is 
heading    NW). 

The  APAIR  MOD  2.6  [Ref.  13:  p.  83]  simulation  uses  a  MAD 
detection  model  that  accounts  for  the  change  in  a 
submarine's  magnetic  moment  (which  is  dependent  on  changes 
in  submarine  heading)  by  using  a  parameter  labeled  DFACTR 
(degradation  factor  for  heading).  In  the  model,  D  (a 
modified  slant  range  at  CPA)  determines  the  probability  of 
detection.  The  value  of  D  is  determined  using  the  following 
relation: 

D    =    DC    (1     -    DFACTR    x    A)  ,  eqn .     4.1 

where  DC  is  the  slant  range  at  CPA  and  A  is  the  acute  angle 
in  decimal  degrees  between  the  submarine  heading  and  an 
East-West  bearing.  The  probability  of  MAD  detection  is 
determined  from  a  table  of  probability  of  detection  against 
slant  range.  A  uniform  (0,  1)  random  number  is  drawn  to 
determine  whether  or  not  the  submarine  is  detected.  The 
average  slant  detection  ranges  (computed  from  Table  IV-4  , 
where  the  vertical  separation  is  200  meters)  for  submarine 
headings  of  North  and  East  are  741  meters  and  545  meters 
respectively.  These  ranges  yield  a  value  of  .003  for 
DFACTR.  The  average  slant  detection  range  from  Table  IV-4 
for  a  submarine  heading  of  NE  is  682  meters;  however,  the 
slant    range    determined    by    a    modified    slant    range    of     545 
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meters  and  a  DFACTR  =  .003  is  643  meters.  If  sin  A  instead 
of  A  is  used  in  Equation  4.1,  then  DFACTR  is  .265  and  the 
slant  detection  for  a  submarine  heading  NE  is  670  meters. 
Since  this  is  only  a  single  data  point  and  there  is  no 
supporting  operational  data,  the  modification  is  not 
proposed  as  one  that  should  be  adopted.  However,  this 
cursory  analysis  does  indicate  a  way  in  which  the  programs 
presented    in    this    thesis    might    be    used   by    others. 

Table  IV-5  lists  lateral  ranges  for  P(det)  equal  to  50 
percent  for  3  submarine/aircraft  heading  combinations.  The 
cross-correlation  and  OPTEVFOR  detection  model  results  show 
the  same  relationship  as  the  results  of  the  square  law 
detection    model. 


Table    IV-5.      Lateral    Ranges    for    P(det)    =   .50    in   Meters    for 
the   Three   Detection   Models. 

Submarine                                    45  0  90 

Aircraft                                    315  90  0 

Cross-Correlation             885  934  754 

Square    Law                               685  730  498 

OPTEVFOR                                     318  358  230 


For  the  detection  ranges  reported  by  OPTEVFOR  [Ref.  5: 
p.  5.1],  the  effect  of  different  headings  was  averaged  out. 
That  is,  measurements  were  taken  from  the  16  possible 
combinations  of  the  4  cardinal  submarine  and  aircraft 
headings  in  equal  numbers  and  then  averaged  to  yield  an 
average  slant  detection  range.  But,  as  shown  in  Tables  IV-4 
and    IV-5,     the    models    show   significant    variability  in    lateral 
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detection  range  for  different   submarine  and/or  aircraft 
headings . 

Figure  4.9  is  included  to  show  the  lateral  range  curves 
when  the  submarine  is  headed  East  and  the  aircraft  is  headed 
North.  These  lateral  range  curves  give  the  minimum  lateral 
detection  ranges  for  the  different  heading  combinations. 
Also,  for  the  cross-correlation  and  square  law  detection 
models,  the  lateral  range  curves  are  fairly  symmetric. 


Figure  4.9.  Lateral  Range  Curves  for  the  Submarine  Headed 
East  and  the  Aircraft  Headed  North. 
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D.   SUBMARINE  DISPLACEMENT 

The  submarine  magnetic  dipole  moment  program  within  the 
main  program  is  used  to  calculate  a  submarine's  induced 
magnetic  moments.  The  program  is  based  on  a  model  described 
by  Forrest  [Ref.  9:  pp.  35-38].  The  model  requires 
submarine  displacement  as  an  input.  Table  IV-6  displays 
results  when  the  submarine  displacement  is  doubled  in  each 
succeeding    entry. 


Table    IV-6. Slant    Detection     Ranges     in    Meters    for 
Different    Submarine    Tonnages. 

Displacement        Signal  Slant        Detection   Ranges    in   Meters 

in    tons                Magnitude  Cross-                  Square                  OPTEVFOR 

in   oersted  Correlation        Law 
cm  ^ 

236 
297 
376 
472 
597 
753 


As  can  be  seen  from  column  two  in  Table  IV-6,  the  dipole 
moment  is  proportional  to  the  displacement.  Since  the  three 
detection  models  give  a  slant  detection  range  that  is 
proportional   to    the  cube    root  of    the   dipole    moment,    doubling 

the     submarine     displacement      should     multiply     the     slant 

1  /  3 
detection    range    by     2     ' J    (1.26).     This     is     confirmed     by 

comparing    the    slant   detection    ranges    between    the    entries    in 

Table    IV-6.       Doubling    the   displacement    multiplies    the    slant 

detection    range    by    1.24    for    the    cross-correlation   and    square 


1000 
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4000 
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8000 
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law    detection    models    and,     as    expected,     by       1.26    for    the 

OPTEVFOR   detection    model. 

Table    IV-7    lists    the    displacement    in    tons  of    selected 

Soviet   submarines.      The   values    were    taken  from   Combat    Fleets 

of    the    World    1982/1983    [Ref.    14:    pp.    602-614].      This    table 

Table  IV-7.   Selected  Soviet  Submarine  Displacements. 

Class  Displacement  in  Tons 

Typhoon  25-30,000 

Delta  III  10,500-13,250 

Yankee  8,000-9,600 

Echo  II  5,000-6,000 

Victor  I  4,300-5,100 

Charlie  I  4,000-4,900 

Tango  3,000-3,700 

Foxtrot  1,950-2,400 

Whiskey  1 ,080-1 ,450 

is  presented  solely  for  the  purpose  of  the  information  it 
contains.  The  submarine  magnetic  dipole  moment  program 
should  not  be  expected  to  give  accurate  estimates  of  these 
submarine's  induced  magnetic  moments,  since  the  program  uses 
a  value  that  relates  displacement  to  magnetic  moment  that  is 
based  on    submarines    of    World   War    II. 

E.   VERTICAL  SEPARATION 

Figure  4.10  shows  three  lateral  range  curves  for  a 
vertical  separation  of  500  meters.  The  OPTEVFOR  detection 
model  lateral  range  curve  shows  only  a  slight  detection 
probability  even  when  the  aircraft  passes  directly  over  the 
submarine.  The  cross-correlation  and  square  law  detection 
model    lateral    range    curves       show    an    increase     in    lateral 
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detection  range  over  the  base  case.  The  dip  in  the  lateral 
range  curves,  for  each  of  these  models,  suggests  the  complex 
variation  of  the  magnetic  signal  with  lateral  range. 


o 


-1000 


METERS 


Figure    4.10.      Lateral    Range    Curves    for    a   Vertical    Separation 
of    500    Meters. 


Table  IV-8  lists  the  lateral  detection  ranges  for 
different  vertical  separations.  It  should  be  kept  in  mind 
that  these  values  are  for  a  single  geographic  location; 
consequently,  they  may  not  be  representative  of  other 
locations.     Note    that    both    the    cross-correlation    and    square 
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re  Law 

OPTEVFOR 

614 

360 

685 

31  8 

720 

225 

724 

-- 

699 

629 

262 

law  detection  models  lateral  detection  ranges  increase  with 
an  increase  in  vertical  separation  until  about  500  meters. 

Table  IV-8.   Lateral  Detection  Ranges  for  Different  Vertical 
Separations. 

Vertical  Separation     Lateral  Detection  Range  in  Meters 
in  meters  Cross-      Sq 

Correlation 

100  804 

200  885 

300  942 

400  974 

500  980 

600  974 

700  936 

No  longer  attains  a  probability  of  detection  equal  to  50 
percent. 

A  factor  related  to  vertical  separation  is  the  effect  of 

ocean  wave  noise  on  a  MAD  system.   As  the  altitude  of  a 

magnetometer  is  decreased,  the  magnitude  of  the  ocean  wave 

noise  increases.   Because  of  the  rate  of  this  increase,  for 

a  given  submarine  and  submarine  depth  there  is  a  minimum 

altitude  at  which  an  aircraft  should  prosecute  a  submarine 

using  MAD.  Further  investigation  using  an  ocean  wave  noise 

model  might  be  valuable. 
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V.   CONCLUSIONS 

This  thesis  has  presented  a  comparison  of  three  MAD 
detection  models.  The  cross -corre la tion  detection  model, 
which  models  an  optimum  detector  under  the  conditions  of  the 
detection  model,  yields  the  maximum  detection  range  for  a 
set  of  given  conditions.  The  square  law  detection  model 
does  not  describe  an  optimum  detector  under  the  conditions 
of  the  model  and  yields  shorter  detection  ranges.  In  the 
stochastic  mode,  with  an  appropriate  choice  for  the 
parameter  AL  that  determines  the  standard  deviation,  the 
lateral  range  curves  for  the  OPTEVFOR  detection  model  become 
similar  to  the  other  two  detection  models.  Detection  ranges 
for  the  OPTEVFOR  detection  model  depend  on  the  choice  for 
the  Operator  Recognition  Factor  (ORF).  With  a  value  of  3 
for  the  ORF,  it  yields  the  shortest  detection  ranges. 
Adjusting  the  magnetic  noise  level  by  an  amount  proportional 
to  the  effective  ORF,  the  cross-correlation  and  square  law 
models  can  be  used  to  describe  the  performance  of  an 
operator    or    an   automatic    detection    system. 

The  magnetic  signal  and  magnetic  signal  plus  noise 
traces  appear  to  have  the  characteristics  of  actual  signal 
traces.  This  suggests  that  the  signal  and  noise  models, 
which  are  the  basis  for  the  cross-correlation  and  square  law 
detection   models,    might    be   useful    for    training    purposes. 
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Variations  on  a  set  of  base  case  parameters  were  used  to 
show  relative  changes  in  the  detection  models.  The 
parameters  included:  magnetic  noise,  submarine  and  aircraft 
magnetic  headings,  submarine  displacement,  and  vertical 
separation.  Significant  results  were  the  large  asymmetry  of 
the  lateral  range  curves  under  certain  conditions  and  the 
variation  of  the  magnetic  signal  as  shown  by  the  changes  in 
vertical    separation. 

The  FORTRAN  and  BASIC  programs,  along  with  an  input 
parameter  discussion,  are  included  to  facilitate  the  use  of 
the  three  MAD  detection  models  as  they  are  implemented  by 
the    programs. 


42 


APPENDIX    A 


LU 


o 


cnjol 

•  v> 

u_ 

a. 

U-Q 

001 

ooO 

*•■»«. 

!-•»- 

*m. 

1- 

LULL 

XOOD 

u_  » 

<• 

0 

"-•> 

XjSOC 

•>• 

• 

Q. 

l-UJ 

9»  9»  m 

— >    * *• 

r» 

Z. 

*-•*- 

-«og-« 

11  Q_ 

u_ 

►-< 

-JQ- 

O^co 

— » 

OO 

•h 

»-*0 

»  •»  » 

» 

z— 

r*» 

O 

CD 

r*^*M 

• 

<Q 

»T 

< 

<io 

I—-*    •» 

UJ 

Q. 

I— 

2. 

cuz 

»    *fw 

21 

h-» 

•— • 

o< 

-*C\J> 

»— < 

QC    - 

•rm 

<t 

QC 

l-O    * 

K 

LUt- 

•  • 

a.    » 

•»  «^y 

->C\| 

r-O 

■«. 

2 

— i— ■ «> 

h- 

1— 

u_—« 

a 

U-<1 

U.l_)    •> 

Z 

•         K 

•■u. 

cc 

O  — 

•  *> 

H* 

u. 

-CUJ 

oox 

as. 

►-J 

h-t_> 

V) 

UoC 

o   *  •*» 

u 

voco 

•UJ 

C3C 

»-l<t 

incvj^o 

•»»» 

niuj 

UJ 

l-O 

— >o.o 

O 

00 

•a. 

H 

<or 

o  •«  *io 

LU 

<-<a. 

vOOO 

UJ 

t-co 

XHN* 

QC 

U_. 

U_ 

3! 

3 

->xz 

a. 

•»  •> 

•»-J 

< 

Q.    •> 

—    •>  »»h- 

-  o> 

p-UJ 

cc 

zz 

OO-rvj  » 

X 

II  -* 

rsjet 

< 

oo 

o  -x— 

< 

-•01— 

1— 

a. 

o»-i 

mm  »o 

2: 

r-z  » 

»* 

i- 

'•QhO 

•<• 

m  + 

LU 

UJ< 

it  -xm 

q: 

in     — 

••N. 

X 

a. -j 

•CO    ■«■' 

LU 

-<ct:o 

>«       » 

>- 

»-UJ 

-»  •coin 

:c 

ILQU 

u_ro 

LU 

oc 

o<xa 

O 

*X  — 

•»  • 

00 

X 

acoc 

o  *  •>  » 

_i 

-       O 

o>o 

H 

<-H 

oo 

inOcr»_j 

•  a. 

— •— t 

3 

H- 

U.O 

-J 

—  »x< 

1 

it    •»» 

l-U. 

Q- 

• 

—1 

CNiCC    »  •» 

>•  • 

•>  » 

Z 

-J 

3t00 

•» 

O   •O^U. 

CC 

Q    •>—* 

m-  — 

>— t 

<5 

<00 

f^ 

-X.2.QC 

a 

«Ji-»f-» 

•UJau 

^ 

au 

m 

— .    •»  •<_» 

QC 

UJ    «H- 

HrfVI    • 

i. 

cc 

OCX. 

UJ 

o-»ro  » 

a: 

wO    *— 

»a.£x.u-\ 

< 

UJ 

uu 

•» 

O    »-S£Z 

UJ 

U.CVJ-    • 

•3-* 

QC 

>— 

ex 

au 

mo  »  •• 

u_o— 

'HUlt 

O 

z 

CLLU 

Z 

—  OCsl«— 

»- 

X   »ZwlU   ■> 

O 

»-• 

X 

Z 

^O^Q. 

3 

H--  QC* 

.  (_         -    *» 

QC 

CCH 

» 

Q    •»   »   » 

Q. 

OC  II       Z    *_l»-iO0 

a. 

^    0 

LU         • 

< 

»Hf> ll/» 

—z 

<tO,QCH-'-*UJX    • 

ir>     ^ 

►-CCCO 

< 

—x*:-j 

«JO>— 1 

UJ       \—  m 

•  cccx-in 

UJ 

U>     < 

JUJ 

•» 

QJ     •»    •»    •• 

•• 

m    O-J    "!*••    •    <-• 

j_ 

i—a.         —• 

U.U.UJ 

2. 

U^>m;w 

U)    •. 

•*«l»     VWU.      »     MU 

*- 

i»uy     «*• 

x    o 

m 

o  »a.z 

1— »•  — 

•v  3E   *>T   ■«•  •     • 

mx»uj     -• 

ozo 

X 

«-<o  •*  •• 

U— ivtO      -4H 

u. 

cr>OQC     0 

ooz 

w  •m»f 

— •  ■->•  .  j—  «(_■  .  f_ 

O 

rg  »u.     <-• 

I-* 

cc 

o-**a 

•0           » 

<»-z 

UJ 

•>  •»  » 

»-h-H-H-|-»— 

•l-H-»-H- 

t- 

r-OX     >t 

oo 

o 

<m>j— ■* 

<<*t<<  <*-»<!<  <t< 

QC 

10  IK     "— 

OOLU« 

UJ 

-i-lSLZ 

zrssrz- 

•ZZZ2 

< 

11  ox    a 

t-tl-l- 

►- 

<   *  .  • 

acacQccZacacoocQCacct. 

a. 

QCO>           *»• 

luo 

z 

UJi    li    If 

uuuouu* 

uuuu 

U«— »*—      uj 

</)(_}  UJ 

•"* 

ac— ujz 

U_U.u-U.U-U_    •U-U-U.U. 

\rt 

-J     a. 

•-"»     t— 

*** 

* 

*—> 

a. 

J.M-UJ 

J> 

z 

KUQ 

0 

h- 

hH 

000000 

0000 

^-ivioi^in^j 

r><juu><-4 

— •  tJ 
-•  o 
x    m 


o 
# 


o 
o- 


o 


QC 

o 
o 

X 

QC 
LU 

a. 

CO 

.< 


LU 

h- 

oa 
•  a. 

00 
00 

vU 

moc 


o 
in 


u_»^  »OHX        • 

o-f\xou.  +  j 
•— r»*^_j'-«X   • 

^QhlLlL+f  Z 

II  <l  ||  «  ||  XCNi«- 
<— lUJf-l  II  U.  M    II  U» 

1-aoxiJ.ZM 


0000 


000 


00 


win    u 

UJ  LU_JLU 
HUOULW 
•-•l-Z*  J 

30H-I-U. 
Z  II 

ur-t— 

Uh3 
Q- 
Z 

O      "-• 

o 
o 


I— LU 

— o 
o^ 

-^X  I  Cv— 

o— x»->o 
^o-?!-o<r 

QULU.XO 
*l  II  II      < 

ujx'^^~uj 

axaua 

a. 
z 


43 


OUJ 

— I -J 

z 
o< 

oa 

OCD< 
•<f 

OS 

u 

•Q 
<< 
<5I 


ULi 

-J 

e> 

z 
«t 

Ol 
I— I 
Q 


U 

1-1        e? 


<       a 

x       h- 
a       < 

-^  U.       UL       -»^ 

O  Q.tfQ«lOO 

»    mui     soct  ••  • 


-3       —  — »-?—» 


•»    .21/1  •»«— 
--*     '-<0«-'Z 


Q 
iJJ 
UJ 

a. 

CO 

Q 

*£ 


LU 

to 
a: 
_}-« 

UICJ 

o 


a 

LU 
UJ 

a. 
<-o 

u 


t-Hi— 10 

003 


LU 
LU 

a. 
00 


to 


a_ 

a 

o 

> 
a 

LU 

<I 


X 

o 

a. 

Q. 


a. 

2E 


QNO|-aiN-40QQSSMOWlJ  I  XX 

<xu.i-z>-  11  11     <<*-jaco  11  11  -)->  11  •» 

uj  II  cjCJ-J»-,*>Jf— ujuj  it   11   II  fii-i**  11  -^i-i 

OO  —1 


«-- H(M  »uj_i  »a: 
cicouj  *-•«_>  a>r«-j 

oc    uoa-c«-zo.«s 

OO       <LU»-it-i       —4 


-JX  I    I  -J<mZ      <I  II 

«l  11  II  11  «■*•  11  utcjr-uji~i 

iwiu-a.c?'-'u_jt^_>a:o 


LU 
CO 

a: 


o    o 


LU 

-I 


a: 

LU 

> 

LU 

-J 

CO 

as 
o«a 

CM 

mat 

O 
OX 


«M  .diCOO 

CJ  HQU       •-»       Qi 

v/#CV   •  CJ 

**r\j     00 

HHJ03    •!— 

>;>  •.*      u>z     *-  » 

I     I  »-*    »  »^UJ       OQ. 

rD(M    rsj— "—  jsoujfnz:<o 

CJCJ       ^»  04  »^-»VJV>HJ  (_!«••       CUU. 

0««*I-03>T     0030 
^Zt/XOOO^J-LUrn    »H-«^ 
•O    •0C   »mQI"   •O.H-J   •'Ohh    » 
xrcO«TO>TClCVLJ«0      h- C3>TUjz»r 

LU       <\|       C\J(\|       LULUO— 1       << 

Qy)QOQ»-JI-fl£3QQ<rQ 
<t      «4I|**II   II  — l»-i       ||       «i—      <l 

ijjf—  uj^uji— K\i«i|XVHMI-  uju.h-  <JJ 
QC3tXU0L-*^OJE»-'a3CX,»-'3IX 

u-  a.        a. 

*~4  \^J  •-*  ^* 


cv      CJ 


44 


LU 

I- 
< 

•■* 

o 

cc 
a. 

a. 
<i 


o 


oo         nO 

— .— i        m 

<T*T  CJ_> 

— QCtthZ 
QUO« 
QQSSOI- 
<«CO  Z 
ujuj  ||    ||         |_J 

oicr.<icu     o 


CC 

o     O^ 

h-LU 
Q      O  — 

J«lOO 

LU«*        <-« 
•»»«■•«-»    •« 

u_oo*^ 
m  •*- 

x  *o 
h-sruji_i 

oc— »  •< 

<      <UJ 

h-UJU. 

a. 
Z 


O 


CO 

o 
o 


CO 

O 

# 
o 


LU 

I- 
< 

»-< 
X 

o 
a: 
a. 
a. 

< 

CL 
□ 

co 


LU 


O 

cn 
m 


CO 


CO 


CO 

O 
O 


LJ 

3£ 


<oooo 

^•4.         t-t^i-^ 

<£—- >»  »  •  • 

<  »o  UJ 


m      •     » 

LUOLUsO 
LJ— >0_>*- 

hZOZ     OOOGCw  «<i<<zm 

Qi- tO*^CU    •    •    *LU       <LUUJCLJi—    • 

Ot—  M-HOOOQ.Q<  cCcC  OC"t—  I*" 
Z  It  Z»-»  •»  H   M       <*-  Z  II 

LJf-tCJOL^«>,*Jh-CUU.  CJ»— • 

u. 

u     o  *■«  cj 

n»     >T  in 

Ml       HI  u> 


uj     «■• 

<      — 

-I     z 

COf-ICO 

oz 

LU^-I* 

-j  -m 

•   •  >-»     t-« 

hhqqui 

II    II        <5  l| 

t^mr—  ujct> 
a. 


CO 

O 

u 

* 


o 

z 

•— « 
CO 

* 


co-*~*Z 
■R-cozcm 

ttl^Vj  — 

cozcoo 

a.cjviw—2. 
♦****+ 

sujzrujs: 
ii  ii  ii  ii  ii  ii 


a: 

LU 

CD 

2. 

<M 

3 

* 

Z 

* 

"~* 

CO 

x 

LU 

X 

O 

^ 

z 

o 

< 

• 

QC 

o 


II  -l_iZ«-. 
CM'M.^tLJCt 


o 

lit 


H-.U 


l_)*->OCJ 


45 


< 

Z  —  LL 

•-4C0Q 

|b 

LULULU 

* 

0 

COJ.LJLU 

■*c 

0 

— 

Ot-«IX 

H- 

•• 

Z 

Oil— 

■• 

fM 

1- 

CCOOH- 

-5 

< 

•k 

U^       <£ 

•> 

■k 

LU 

LL.CT  <!•-• 

—1 

r- 1 

•» 

3 

•» 

i«i 

—1 

-J»-Q— 

Z 

•» 

-J 

— 

<aujii<i 

«k 

1— < 

•» 

»-* 

zac<i2. 

CD 

z 

.-4 

•— 

(^        0 

4^ 

•» 

CO 

-c 

•— "OOf-" 1 

•■ 

^H 

•I 

h- 

OOOOLJCO 

1— 1 

-> 

00 

•» 

—II— 

ii 

•k 

•^ 

z 

-» 

LU<L_>  II 

•> 

t-4 

CO 

•» 

— 

1-^ 

-L       LU 

CD 

CD 

•^ 

— 

t— 4 

LU 

— 

h-l—  >:> 

O 

•> 

— z 

LL 

*z 

«» 

O 

1-1      LU 

•k 

O 

CO  1-4 

»—* 

•ft 

(N 

•i 

LL      OQ 

0 

Z 

—•00 

— z 

O 

Q 

«-» 

O      O 

CO 

•» 

z# 

U_H-I 

•ft 

•a 

*-4 

•    a 

•» 

O 

1    1  — 

—  CO 

o 

-4 

•M 

ujooLua: 

m 

-5 

O0Q 

z«- 

00 

— 

a 

LO 

Q_IX< 

Q 

■M 

■S-  — 

*-«— « 

•» 

rxi 

m 

Q 

Z3LUJ-Q 

•» 

O 

—CO 

coa 

m 

•» 

r* 

•» 

H-Q      Z 

2: 

CDZ 

QO 

^J-  — 

a 

LU 

;> 

*■ » 

mQZ< 

•ft 

•>  » 

—O 

—CO 

•ft 

•i 

m 

<— • 

Z^1-1*- 

CL 

moo 

Z  1 

ao 

3: 

m 

2: 

•m 

O            CO 

•» 

Q2 

n^'~ 

—0 

■• 

Q 

«» 

LU 

<N 

OLU 

LL 

«.  •> 

co— 

z  1 

a. 

■i 

— < 

—1 

Q 

Za:oi 

■h 

Q-  — 

1  < 

1—4.^ 

m 

0* 

00 

1-^ 

» 

OZK 

< 

►O 

—  1 

CO  — 

LL 

J 

•« 

LU 

— 

'ilLU<i— 

•» 

u_o 

<o 

1  0 

-^ 

•» 

• 

.— < 

h* 

iiajE 

0 

•^n 

1  0 

—0 

z 

<x 

-J 

^ 

h- 

— 

H-l-      — 

0 

<— 

0— ' 

0— 

* 

•» 

<l 

«» 

3 

—4 

—1 

4» 

*z 

—    oz 

oz 

0 

o 

» 

<i 

u. 

a 

(/>^4C|,UJ 

LU 

o»— 

O        w-* 

«*»-^ 

«£ 

o 

LU 

> 

1— 

«• 

LULJ<XV> 

m 

0  » 

a.     zco 

ZCO 

1 

•» 

C£ 

•ft 

_> 

*» 

^-^-H^JJ•-* 

ISI 

••— • 

»—*-«» 

t-..JJ. 

—4 

oj 

(-1 

^ 

UI 

«— < 

JKt-U 

•> 

OJO 

CJ<LCO— 

to— 

-> 

•» 

» 

m 

_ 

CLU<Z 

O 

-1  -O 

•  1  4S-CQ- 

•*u» 

* 

rsl 

z 

0 

0 

0 

SLU_J 

X 

•o»-« 

1^0—— O— mo 

m 

•» 

00 

< 

?»; 

Oh-     2: 

■w 

N^X  — 

•OCDCOOLUCO*  -J 

~*cocr» 

a. 

— 

2. 

00— 

OLU— O 

2: 

•»   m(J) 

0—0- 

— 0* 

1 

u. 

"f"*_j 

— 

f-1 

4^0 

CJ^iLJ 

0 

-»oo 

XC)C0CJOC0V->O-H 

i— iwr* 

a.—*—* 

— ^« 

m 

LU 

U 

*if 

LU       — ^ 

*/) 

•  ••  » 

"— uu*  l-icj*  a. 

LU 

rMl/KT.^, 

lj^-j 

<JU    •iU'U'J'U 

r-  — » 

»  •» 

4C-J       *4 

«INIJ 

UJtJV^— «««»«-»—  > 

'# 

— ■  -».^»-» 

.^.^-« 

Z^Z3_JLJQtt 

o^r* 

t-KlOQC 

LU 

*>   »   « 

l-*«.0*#Q 

0 

oooo 

ooo    r- 

—  II  CO  ||  X  +  Q- 

00>O  11  — 

I—  ZLU 

z 

soo 

i->t  |i  ii-i 

-4^rnr*-z 

XLU       - 

.       (^       LUCQ    • 

H-tLULU 

_00»-00 

•— 1 

X   • 

i—  CDC2U.QZ    • 

* 

*•»•»■• 

•»  ••  •*:* 

—        _JIAJ_J_J-J-JLJP- 

H--3 

0>-4<J^     • 

K- 

(X      C> 

LJ*^'■^,■^"," 

»— I-4VJ 

1 

•4*  «*•  ^  -4-  »-•  «r  -4-  -r  -S.O  oo  u.  o  -j  * 

--HI-IZQC  — 

o«z 

oc^OLU—io: 

3 

LU        • 

ac  co  co  co  coco  co  «^o 

«Mww(^w — wco  «a 

:»-»h-<co  <o<  *-*o- LUGOac1-" 

03      -JQ.O 

O 

O      >-i 

OO—OO— V 

.    • 

QQOQ* 

QQQiiM  | 

II  lf>OXU-JOI-      1—  LT\3>-a. 

OOLU      K- 

0C 

LU^i-) 

— IOO  1  UU  1  CLCN 

•<t«l<  II    II    II  00 

ZH'h 

ZOO 

V53WJU 

03 

H— «    - 

-J  II    II    II   II 

II  II  II 

II 

UJLUbUUJ    II 

aituaiTCOV 

•^u 

UJU.LJ      v_ir— ^. 

h-      <IUJ 

_> 

4CLUi 

«*VJOCJ  «-»!-*  «-»•—■ 

rvj 

aasK^ixaccQ^uzu 

OQ--SO     uviai 

V)      <z> 

CO 

•-•eta: 

CJCQ  "J  ^  CO  "J  <£  ^  «* 

»- 

t-4LU       O 

* 

-j 

-LJ-LU— -£. 

OLJ 

0 

f-t—  OCOI— 

r» 

au 

Ui 

ui 

uuuuuuu    o 


46 


< 

•— i 

a 


UJ 

X  fH 

o 

u.  —  •— 

UJ  o 

z  •  o 

U  fN  O 

mm  Q  l-l 

h-  »  w 

O  *-«  CNJ 

UJ  Q  Q 

—                                              uj  r»  — 

o                                            a  >  o 

— •                     5£                                         U,  S  O 

»               ««•                              o  «*  — • 

c\i                  o  ^  •*. 

a.                C9                               co  </>  -^ 

+                           a  *  a 

Z               — »                              ac  ■— t  »on 

h-              m                            o_  x  r-> 

in                       ••              2:  •»  >  •» 

•                     r»               ■+•                            uj  -o  »m 

r\JO                     no     ^c         x                                 x  >  —ico 

**                     m     #         —                                »—  »  x  •• 

*•—<_>     o    -*«r     —        4&                                     •  x  »o 

-*ocr     co     men     x.        >—                               coco  »  -»_j 

en     o<                  •  rvi     w         ||                                z_j  o  o  •>               — 

ou     ir+o     uj     o-^»v^.        —                               acuj  co  om                ** 

•#     ara>-    h-    #<— s:»-    rm                           rjQ  «■»  o<             co 

1 0         -to     o     ict—<— •     — •+                                 ujz.  ca  —  2:             — 1-* 

—    •      0<                     "—  C  II  XUJ  II  XUJ                                     Ct  O  XZ.                COCO 

1  »— uj*^    a."-  >»    — ««j  -j  -x.  a.«M—»  _j     ^*c     ^^t        o               cute  u»  »ovu           iv> 

•h-a.     *-••»    ?*    §«•<—     01—  i-iuf- »-«uli     -*                  ^u»  cu.^.j»     <-»     ac<_>— » 

ZO      3!   •—  yQ      «^<Z      >saruJC?UJ3:,-,— ItTV      H- *0      K3      O                   '-•UJ  UJ  »  •  •     o     *cos: 

•-.Z      »-«r\_j~-      IIOIO      »Z  HZ  II  II -J          Z          Z2      +                    H-X  Z  STOCvJ     CO     -H«  «~ 

«#      H-H-   «UJuj»-«-J*-»XO<\l    •— •   «Nm<0      OO      Oi-iUJOUJ               3h  •-»  COl£      <—     ^-*»-o 

«i«-ujh-^>    u.ihn  »<n»-"»-<a:a.tjt-i    <_>«->    u^jwj           cj  h-  cc           f— c\j*it<t.j 

OOOQt-ZHZO  II  **Wh— «-                                      zz  II  zz         etz  =>  UJ             QiX>0*O-f 

•  •  ».j<x  «a:»-»o»—  11  11  11  nu.ou_                               OMO«a        coo  o  o           omi>^_js: 

<nf\|OU_3:"9  3»— <«WOOOmOm                                            OI-CO|-Ij           3»-«  Qi  LU»J_J      CO—  |  ^U.Z 

II   II   II   HtXw      Z                                                                                            Z      Zl-Q      COH-  CO  H-<k       HUJII   ||    M   II 

-njwiuu.     i_j<_>                                                                    cj     <-jujz        o  _?  zujlu     (\j-oooiPi 

^<ivixu.i-(     00                                                                          o     <w»a.uj     touj  to  "-^cxqc     ac^^^a:**. 

— •>— 

J.UJ 

OO  CJ         000     o 

04         ni  u»         'uwv     1-* 


-"» 

>T 

* 

►H 

0 

fc 

—  ~3 

• 

»-•  1 

Z^ 

U. 

*cu 

LJW 

0* 

fO 

COO 

UJ 

"♦■CO 

000-5 

■-•  1 

•<»-« 

CD^-4 

^aoo 

•*<■ 

1  r-U»»N 

CJ* 

— »•          M 

47 


X> 
CO 

CL 
O 

h- 
O 

UJ 


O 
UL 

> 


CQ 
< 
CQ 
i-l 

cc 

Q. 


Z 


3D 
□ 

Q- 


O 

— ■ 

I 

CO 
CQ 

< 
© 


cc 
or 
co 


m 
«~co 
0>* 

+  > 

a. 
-»cj 

">.  • 

CPCVJ 

+  t- 

IX 

oo 
•co 

"■«  I 
II  II 

mo> 
coj» 


CO 


C0>-CO?- 

O  II  O  II 

or—cc— 

U.LUU.UJ 

•— fw  —  — ^ 
-•     ogee 
JQ-IQD 
-J      _J      h-Q 

<t       <%       LU^. 

CJ         <w>  U.LJ 


(J 


z 

O 

< 

z 
< 

UJ 

-J 

z    — 

<    *: 


t—    .> 

•i 

CO     -J 
z    — 

(X 

3     UJ 


UJ  UL^» 

z  » 

1-4  LUX> 

I-  z  * 

CO  o 

-D  <X_J 

00  CQ<1 

co  cot 


ac 

— » 

CC 

u 

«■* 

a 

^ 

> 

z 

LU 

rn 

UJ 

X 

X 

X 

h- 

+ 

CM 

h- 

cc 

X 

CC 

a 

••» 

0 

u_ 

* 

> 

u. 

*• 

+ 

Nl 

M 

—4 

*^ 

X 

z 

U. 

* 

> 

Z 

X 

> 

•— • 

UJ 

•» 

♦ 

ID 

> 

m 

i-* 

X 

0 

— » 

<J 

(NJ 

• 

rvi 

r>l 

— 

X 

> 

M 

U. 

CJ 

J- 

•» 

I  *"* 

"V 

<—• 

»4 

00 

»• 

— * 

*— 

~ 

00 

>- 

Z 

X. 

X 

o> 

— • 

^: 

•» 

o 

CC 

«-• 

•» 

•  >» 

rgj" 

ac 

^ 

1— 

3 

•-4 

eg 

*-*> 

O  I 

r> 

<— ' 

a 

H- 

2! 

O 

II  "N, 

*  II 

>— 

co 

—o— 

LU 

CC 

•» 

>■ 

>->- 

UJ 

a 

^    > 

(X 

a 

1-4 

0 

+ 

CC 

CC 

*— •  •» 

z 

<jj 

— .    • 

rM<^ 

a. 

-*    _j 

UJ 

^■» 

«» 

^H 

Ul 

LU 

"-• 

+  <-»— 

■£. 

cv 

u  t«-»r»o  lajvju  vj'uj—*  ut 

*c 

_J    »C\1C 

) 

«•-« 

UJ 

C3 

•  Or-nnrgaj^OJfr    • 

— 

UJ 

*0-£C\ 

1 

H- 

Z 

m 

ooinoarnr«'(Nim>o 

H 

Z 

ID   •< 

UJ 

UJ 

3 

M 

> 

•  _jmc\JorsjO*»-*+    • 

3 

►— t 

— CJH-LJ-J 

-J 

Ul 

H- 

h-*.HU<-tmWOUK 

U 

1- 

h- LU<»— Z 

■z^. 

cc 

3 

OH*ir>coo«*" 

-*oooz 

Qd 

3 

a:  •  n  a 

i«-iO»-»oC 

CO 

O 

•  C£      •     •     •     • 

•    »^«    •  QC 

CO 

O 

&^-iU>t- 

•  1-3 

3 

cC 

-J 

XONOOH 

OO  1  X3 

r> 

X 

U1-— 

ZOZI-Q 

00 

CQ 

< 

X  wCO  II    II   II    II 

II    II  >— H-Q 

00 

CO 

II  u_ 

CJ 

II  UU-*4£ 

-J 

UJ 

11  u_  11  cj^^rgi— 1 

iMm  ti  u»uj^: 

_> 

5t>-i 

W~>CJCCUJ 

OO 

OO 

CC 

j-<-*>-e>e>ox 

jl-a->-— •CI.LU 

00 

CO 

•— • 

M 

— 

J. 

t— 

H- 

o 

CJ 

i-H 

»>* 

uuu 


CJCJCJ 


ouu 


48 


UJ 

Q 
O 

2. 

cC 

o 


Q. 

O 

«•  UJ 

a.  h- 

•» 

O  CC 

•  t_l 


CC  UJ  ■> 

to  » 

^s  U_  > 

—  o 


in 

>■  or 

>  ft 

Q-  «* 

Of 
in  <— 

CJ  -»  ft 

*  .   >  3 
■T              »  + 

»v  o     ar>-  > 

PI  t-t  o     *■»*  I 

CjT  <*■  -*       #»~ 

•  vO  3—  O 

«M  ^^  U l        f  CJ>  • 

or         i  m  H->tn4  •  — » 

»         M  <N  Oinonj  it 

— <  >    •       (M      00       0^**^  > 

3  o     cor-r-cr*     cm«-> 

Jfc  U\l<tU,U)V'CJUIt> 

>  •  •*Jumr^rsir-in,TO  I  -*  • 

x         •«— CMn^cvjo  •  •ftQ.     uj    uj  • 
►—  >»»-»PtcouJOi*-oijEx<_J:3     _?•** 
j     m  «r-  »rno  h  h  \uy-~Z.     zoz 
•o  «o  «»^  •  •  •— «^  ti  a«— <o*-«  «cc 
_j        x  «o  I  -»  I  -»>o.>>oi-  »I-X3 
<    x<-»^  ii  n  ii  ii  H  _  zoz—'f-a 

uj      ii  u»  ti  ~r\4cn>rmu.  u  »  uu«iu^ 


o    o 


t- 

•k 

O 

»— « 

LU 

— « 

-J 

•> 

00 

h- • 

in 

(, 

CD 

a 

o 

«* 

* 

•■ 

CO 

o> 

X 

o 

-j 

CC 

a 

» 

» 

a. 

-i 
< 

o 

UJ 

•» 

o 

X 

u. 

in 

»— 

CC 

*■» 

'-J 

m 

in 

IS) 

•» 

(_> 

• 

Z 

z. 

•» 

a. 

o 

(X 

•» 

a* 

# 

# 

3 

s: 

_i 

ft 

ft 

H- 

^»» 

» 

^% 

•—» 

LU 

m 

-j 

-^ 

<\J    — 

CC 

a 

<i 

z 

*    > 

u 

*> 

* 

*      • 

UJ 

i 

u. 

u. 

^liCN 

-i 

u. 

u. 

►-»•«»» 

•H 

UJ 

CJ 

O 

+  </)CU 

H 

^ 

UJ   •» 

o—x    >»a 

■D 

l—t 

2 

•*v»CC  CM -»CC 

O 

t— 

a.  •« 

mac*-*  Q.U.J- 

CC 

Z3 

UJZ 

o^»*-»#  1  '-•II  z 

CO 

a 

o 

nOU<0>I      —CC 

Z3 

cc 

UJ_| 

•    •«— 

II  -JCC-ILUID 

OO 

cu 

\-<3. 

O-*  II  O— »w-J— h-O 

~J 

ZUJ 

II   II  X"-  II   II  «*"»uj*i, 

V) 

SSi 

•-.(X 

uacccoaNuaau 

<->uo 


49 


LTP    FNG 

PO(CC) 

FC(SL) 

PD(OP) 

G(  I) 

TN(I) 

-15C0.C 

0.0  54 

C.016 

0.0 

0.0114 

C111C 

-147CC 

0.058 

C.016 

0.000 

0.0140 

-0.0488 

-144CC 

0.063 

C.C16 

c.ccc 

0.0173 

-C.17C9 

-1410.  C 

0.068 

C.017 

C.000 

0.0212 

0.C289 

-1260. C 

0.0  74 

C.017 

C.000 

0.0259 

-0.0225 

-1350. C 

0.081 

C.017 

coco 

0.0316 

-C.G464 

-1220. C 

0.089 

C.018 

0.000 

0.0383 

-0.0529 

-1290.  C 

0.098 

C.018 

c.cco 

C0462 

-0.0C6C 

-1260.  C 

0.110 

C.019 

C.000 

0.0553 

-C0C92 

-1220. C 

0.  123 

C.020 

0.000 

0.0655 

0.0C26 

-12C0.C 

0.129 

C.C21 

COCO 

0.0766 

0.GC74 

-1170. C 

0,159 

C.022 

0.000 

0.0382 

0.0248 

-11*CC 

0.162 

C.C24 

0.000 

0.0996 

0.0406 

-1110. C 

0.210 

C.026 

coco 

C1100 

-0.C866 

-iceo.c 

0.245 

C.029 

0.000 

0.1181 

-0.0688 

-1C50.C 

0.287 

C.022 

c.cco 

0.1229 

-C1846 

-1C20.C 

0.327 

C.036 

cooo 

0.1232 

-CC197 

-990. C 

0.398 

C.043 

o.coo 

0.1185 

-0.0662 

-96C.C 

0.469 

C.051 

c.occ 

C.1C89 

0.C241 

-92C.C 

0.550 

C.063 

0.000 

0.0951 

0.0156 

-sco.c 

0.640 

C.C80 

cooo 

0.0784 

0.0616 

-670. C 

0.7  23 

C.105 

cooo 

C.0604 

0.0129 

-64CC 

0.823 

C.143 

cooo 

0.0427 

0.0901 

-€10. C 

0.898 

C.199 

c.ccc 

C.0265 

0.0252 

-760. C 

0.9  53 

C.282 

ccoo 

0.0127 

0.1639 

-75C.C 

0.9  64 

C.401 

cooo 

0.0016 

0.0252 

-720. C 

0.996 

C.557 

coco 

-0.0069 

0.0169 

-690. C 

l.OCO 

C.733 

cooo 

-0.0130 

0.0202 

-660.  c 

l.OCO 

C.8  86 

0.000 

-0.0170 

-0.1250 

-620. c 

1.0  00 

C.972 

c.occ 

-0.0195 

-C1C27 

-6CC.C 

1.000 

C.997 

cooo 

-0.0207 

-0.0647 

-57C.C 

l.OCO 

l.COO 

ccci 

-C.0210 

0.1751 

-540. C 

1.000 

1.000 

0.0C4 

-C.0207 

0.2297 

-510. C 

l.OCO 

l.OCO 

0.011 

-0.0199 

0.1687 

-46CC 

l.OCO 

l.OCO 

CC27 

-C0189 

0.1127 

-4SCC 

l.OCO 

1.000 

C.C59 

-0.0178 

0.1272 

-420. C 

l.OCO 

1.000 

0.116 

-0.0166 

0.1432 

-390. C 

l.OCO 

1.000 

0.2C1 

CO 

0.1615 

-260. C 

1.000 

1.000 

0.314 

0.0 

-0.1C59 

-22C.C 

l.OCO 

l.COO 

C.444 

0.0 

0.0659 

-3C0.C 

l.OCO 

1.000 

C578 

CO 

C1731 

-27C.C 

l.OCO 

1.0  00 

C7C0 

0.0 

0.0912 

-240.  C 

l.OCO 

l.COO 

C.799 

CO 

C1C2C 

-210.  C 

l.OCO 

1.000 

0.872 

0.0 

-0.0259 

-16C.C 

l.OCO 

1.000 

C.922 

CO 

-0.1278 

-150. C 

1.000 

1.000 

C.953 

CO 

-0.0262 

-12C.C 

l.OCO 

1.000 

0.971 

0.0 

0.0S56 

-9cc 

l.OCO 

l.COO 

C.981 

CO 

0.1121 

-60. C 

1.0  00 

1.000 

C987 

CO 

-C0CC9 

-20.  C 

l.OCO 

1.000 

C989 

0.0 

-0.0580 

O.C 

l.OCO 

l.OCO 

C.99C 

CO 

O.0C91 

20. C 

l.OCO 

1.000 

C.989 

0.0 

-0.075C 

60. C 

l.OCO 

1.000 

C.987 

0.0 

-0.0543 

90. C 

1.0  00 

l.OCO 

C.981 

CO 

0.01C5 

120.  C 

l.OCO 

1.000 

C971 

0.0 

-O.OCCC 

lfC.C 

l.OCO 

l.COO 

C953 

CO 

0.0766 

180. C 

1.000 

1.000 

0.922 

CO 

-0.0412 

21C.C 

l.OCO 

1.000 

C872 

0.0 

0.0262 

240. C 

l.OCO 

1.000 

C799 

CO 

0.0C57 

270. C 

l.OCO 

1.000 

C700 

0.0 

-0.026C 

2CC.C 

l.OCO 

1.000 

0.578 

CO 

-0.1266 

220. C 

1.0  00 

1.000 

0.444 

0.0 

0.1146 

26C.C 

l.OCO 

1.000 

0.314 

0.0 

-0.1252 

39CC 

l.OCO 

l.COO 

C2C1 

CO 

-0.0248 

420.  C 

1.000 

1.000 

0.116 

0.0 

-0.1144 

45C.C 

l.OCO 

1.000 

CC59 

0.0 

-0.0815 
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4€0.C 

1.000 

1.000 

C.C27 

0.0 

0.126C 

I  10. c 

l.OCO 

l.OCO 

0.011 

0.0 

0.OC04 

540. C 

LOGO 

C.990 

0.C04 

CO 

0.2C16 

570. C 

l.OCO 

C.925 

0.001 

0.0 

0.072C 

6CC.C 

l.OCO 

C.765 

C.OOO 

0.0 

-0.12C5 

620. C 

0.996 

C.561 

C.000 

CO 

0.126C 

£60.  C 

0.980 

C.381 

0.000 

0.0 

-C052c 

6SC.C 

0.938 

C.251 

C.CCC 

CO 

-0.01C6 

720. C 

0.8  62 

C.167 

C.OOO 

CO 

CO 

75C.C 

0.761 

C.115 

C.OOO 

0.0 

0.0 

7£CC 

0.650 

CO  82 

C.CCC 

CO 

CO 

610. C 

0.543 

C.062 

C.OOO 

0.0 

0.0 

640. C 

0.447 

C.048 

C.OOO 

CO 

0.0 

270. C 

0.367 

C.039 

COCO 

0.0 

0.0 

SCC.C 

0.301 

C.033 

0.000 

0.0 

0.0 

920. C 

0.249 

C.C29 

coco 

CO 

O.C 

96C.C 

0.207 

C.026 

coco 

CO 

0.0 

990. C 

0.174 

C.023 

coco 

0.0 

0.0 

1C20.C 

0.148 

C.021 

C.CCC 

CO 

CO 

1C50.C 

0.127 

C.020 

0.000 

0.0 

0.0 

1C6C.C 

0.  110 

C.C19 

C.OOO 

0.0 

0.0 

1110. C 

0.096 

C.018 

C.OOO 

CO 

CO 

114CC 

0.0S5 

C.018 

0.000 

0.0 

0.0 

1170. C 

0.076 

C.017 

COCC 

0.0 

0.0 

12C0.C 

0.068 

C.017 

C.OOO 

0.0 

CO 

122CC 

0.062 

C.016 

0.000 

0.0 

0.0 

1260. C 

0.057 

C.C16 

C.CCC 

CO 

CO 

1290.  C 

0.0  52 

C.016 

0.000 

0.0 

CO 

1220. C 

0.048 

CO  15 

0.000 

CO 

0.0 

1250. C 

0.045 

C.015 

C.OOO 

CO 

CO 

132C.C 

0.042 

C.015 

C.OOO 

0.0 

0.0 

14  ICC 

0.040 

C.015 

COCC 

CO 

CO 

1440.  C 

0.038 

C.015 

coco 

0.0 

0.0 

1470.  C 

0.036 

C.015 

C.OOO 

0.0 

0.0 

15C0.C 

0.034 

C.015 

CO 

CO 

CO 
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10  DIM  G<!200)>  DK100),  D2<!100),K<100),  X0C100) 

15  DIM  D5C100) 

20  DEG 

30  PRINT  "MAX  FREQ"; 

40  INPUT  Fl 

50  PRINT  SPRINT  "MAX  FREQ  =  ";F1 

60  T1=1/F1 

70  PRINT  "INTERVAL  TIME  =»  ";:  INPUT  T7 

30  G=T7/2/Tl 

30  H=INT(G) 

100  H=H+INT<2*<:G-H)  ) 

110  M=2*H+1 

120  IF  M>200  THEN  70 

130  PRINT  "INT  TIME  =  ";T7 

140  T7=T1*M 

150  PRINT  "ADJ  INT  TIME  =  "TT7 

1E0  PRINT  "SAMPLE  SIZE  =  ";M 

170  PRINT  :  PRINT  "F/A  RATE  "', 

130  INPUT  F2 

130  PRINT  : PRINT  "F/A  RATE  =  " "F2 

200  P1=F2*<:M-1)*T1/3S00 

210  PRINT  "PF  =  ";P1 

220  PRINT  : PRINT  "INPUT  DIP  ANGLE  ( 1=YES, 0=NQ) ";: INPUT  A 

230  IF  A=0  THEN  2S0 

240  PRINT  "DIP  ANGLE  PHI  "',  :  INPUT  F 

250  GOTO  420 

2B0  DEG  :L1=7S:L2=100 

270  PRINT  '.PRINT  "LATITUDE  "; 

2E0  INPUT  L 

230  PRINT  : PRINT  "LONGITUDE  "• 

300  INPUT  0 

310  PRINT  :  PRINT  "LAT  =  " 'L:  PRINT  "LGN  =»  M:0 

320  F=SINCG-L2)*C0S<L)  :G=COS<:0-L2)  *COS<:L)  :H=SIN<:L) 

330  U=G:V=H 

340  GOSUB  1300 

350  J=J-(30-Ll):G=K*SIN<:j)  :H=K*COS<:  J) 

3E0  U=G:V=F 

370  GOSUB  1300 

330  F=K:G=0:R=-<:COS<:L1)*SIN<J))  :a=-<.COSCLl)*CaSCJ)) 

330  U=H:V=F 

400  GOSUB  1300 

410  F=ATNC2*<SIN<J)/C0SCJ>)) 

420  PRINT  "PHI  -  " !F 

430  PRINT  .'PRINT  "DIPOLE  COURSE  ";:  INPUT  CI :  PRINT  "DIPOLE  SPEED  ";:  INPUT  VI 

440  PRINT  : PRINT  "SENSOR  COURSE  ";: INPUT  C2: PRINT  "SENSOR  SPEED  ";: INPUT  V2 

450  PRINT  "DIPOLE  COURSE  =  "?C1: PRINT  "DIPOLE  SPEED  ■  ";V1 

4S0  PRINT  "SENSOR  COURSE  -  ";C2:PRINT  "SENSOR  SPEED  =»  "  ?V2 

470  U1=V2*SIN<:C2)-V1*SIN<:C1)  :W2=V2*C0Sf.C2)-Vl*C0S<Cl) 

4S0  U=W1 :V=W2 

430  GOSUB  1200 

500  C0=J:W0=K 

510  D3=W0*T1*4.  S3/9 

520  PRINT  "REL  COURSE  =»  ";C0: PRINT  "REL  SPEED  =  ";W0 

530  PRINT  :PRINT  "INPUT  DIPOLE  MOMENT  Cl-YES. 0=NO) " ;  : INPUT  AA 

540  IF  AA=0  THEN  5S0 

550  PRINT  : PRINT  "MAGNITUDE  P  ";: INPUT  P: PRINT  "HOR  ANGLE  w  ";: INPUT  A 

SE0  PRINT  : PRINT  "VERT  ANGLE  OMEGA  ";: INPUT  9 

573  GOTO  340 

533  PRINT  :  PRINT  "INPUT  EARTH  FIELD  (.  1=YE3,  0=NO)  ";:  INPUT  QA 
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590  IF  AA»a  THEN  S20 

603  PRINT  : PR INT  "EARTH  FIELD  "?: INPUT  El 

Sia  GOTO  S33 

S23  E 1  =73000/  SDR  ( 3*C0S  (.  F )  *COS  t  F>  + 1 ) 

S30  PRINT  :  PRINT  "EARTH  FIELD  =»  ";E1 

640  M4=0:M5=0:M5=0 

630  PRINT  : PRINT  "INPUT  PERM  MOMENTS  <1=YES.0=NO)  "?:  INPUT  AA 

SE0  IF  AA=0  THEN  S90 

S73  PRINT  :  PRINT  "LONG  MOMENT  " ',  :  INPUT  M4: PRINT  "TRAN  MOMENT  ""INPUT  M5 

SS0  PRINT  "VERT  MOMENT  ";:  INPUT  MS 

690  PRINT  :  PRINT  "LONG  MOMENT  =»  ";M4:  PRINT  "TRAN  MOMENT  =  "?M5:  PRINT  "VERT  MOMEN 

T  =  "?MS 

700  K1=7.3:K2=1.S:K3=1.5 

713  PRINT  :  PRINT  "DISPLACEMENT  " ', 

723  INPUT  Nl 

733  PRINT  :  PRINT  "DISPLACEMENT  =>  "TNI 

733  NN1=N1 

740  M9=£1*K3-*N1*SIN<:F1  :M3=M9+MS 

750  M3=El*C0S<.F^*Nl*(Klt«C0S(:Cl)*CCS(Cl)+K2*SIN(:Cl)*SIN«:Cl)) 

753  Ml=M4>t.3INiCl)+M5*CaS(C1.1  :M2=M4*COS<:Cl  )-M5-*SINCCl) 

773  M7=ei»<cascF)>H<:Ki-K2?*Ni*siN(Ci)*cas(:ci) ) 

7S0  Ml=M7+Ml:M2=MS+M2 

790  U=M1 :V=M2 

300  GOSUD  1330 

313  A=J:U=M3:V=K 

323  GOSUB  1300 

333  P=K:B=J 

340  PRINT  :PRINT  "P  =  ";P:PRINT  "w  =»  "?A:PRINT  "OMEGA  =  ";B 

350  X=Pl 

3E3  GOSUB  1323 

370  V6=Y 

330  V7=M*< l-2/3/M+Y*SQRC2/9/M) ) *3 

390  PRINT  -.PRINT  "VERT  SEPARATIGN  "*,:  INPUT  Z 

900  PRINT  : PRINT  "VERT  SEPARATION  -  *',Z 

910  PRINT  :PRINT  "NOISE  "', 

920  INPUT  SI 

925  S2=S1*4 

930  PRINT  :  PRINT  "NOISE  =  "?S1 

940  PRINT  : PR  INT  "MAX  LATERAL  RANGE  " * : INPUT  RS 

933  PRINT  :  PRINT  "NUMBER  OF  INCREMENTS  " ' '■  INPUT  N7 

953  PRINT  -.PRINT  "MAX  LATERAL  RANGE  =  "',RS: PRINT  "NUMBER  OF  INCREMENTS  =»  ";N7 

954  PRINT  : PRINT  "  ORF  ";: INPUT  ORF 
9ES  PRINT  :PRINT  "  ORF  =  ";ORF 

970  D4=«3/N7:NS=2*N7 

974  PRINT  : PRINT  "  ALPHA  "?: INPUT  AL 

975  PRINT  -.PRINT  "  ALPHA  =  "  ?AL 
9S0  L9=-RS 

990  FOR  E=0  TO  N8 
1000  X0=L9:X0CE)=»L9 
1010  GOSUB  1S20 
1020  L3=t-3+D4 
.1030  GOSUB  1300 
1040  NEXT  E 
1044  GOTO  1230 
1050  GRAPHICS  3: COLOR  1 
10S0  XX=INTC310/NS) 
10S5  X0(0)=0 
1070  FOR  1=0  TO  N8 
1080  D1<:I)»INTCC1-O1(I))*1S0) 
1090  D2(I)=INTCC1-D2CI))*1S3) 
1095  D5a)=»INT(:Cl-05CI)  )*153) 

1100  X0<:i+i)=x0<:  n+xx 

1110  NEXT  I 
1123  PLOT  X0<0),D1<:0) 
1130  FOR  1=1  TO  N8 
1140  DRAWTO  X0CI)tD1CI) 
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115(3  NEXT  I 

nea  plot  xaca),  D2C0) 

117(3  FOR  1=1  TO  N8 

USB  DRAUTO  xo<:  n ,  D2<:  I) 

1120  NEXT  I 

1192  plot  xa<:a),D5(0) 

1134  FOR  1=1  TO  N8 

1135  DRAUTO  Xa<:i),D5CI) 
1133  NEXT  I 

1200  PRINT  "PD  FOR  X  FROM  "?-R8?"  TO  " ; RS 

1230  GOTO  1610 

1230  PRINT  "FOR  HARD  COPY  ENTER  ' 1'" ?: INPUT  CC 

1225  IF  CCOl  THEN  GOTO  1050 

1330  LPRINT  "MAX  FREQ  =  "?F1 

1310  LPRINT  "ADJ  INT  TIME  =  " 7T7 

1323  LPRINT  "SAMPLE  SIZE  «  "?M 

1333  LPRINT  "F/A  RATE  =  "?F2 

1335  LPRINT  "PF  =  "?P1 

1340  LPRINT  "LAT  =  ".?L 

1333  LPRINT  "LQN  =  "?Q 

1353  LPRINT  "PHI  =  "?F 

1373  LPRINT  "DIPOLE  COURSE  =  ";C1 

1330  LPRINT  "DIPOLE  SPEED  =  "!V1 

1330  LPRINT  "SENSOR  COURSE  =  " ;C2 

1400  LPRINT  "SENSOR  SPEED  =  "?V2 

1410  LPRINT  "REL  COURSE  =  "?C3 

1423  LPRINT  "REL  SPEED  ="?U0 

1430  LPRINT  "EARTH  FIELD  =  " ?E1 

1440  LPRINT  "LONG  MOMENT  =  ";M4 

1450  LPRINT  "TRAN  MOMENT  =  ";M5 

1460  LPRINT  "VERT  MOMENT  =  ";M5 

1470  LPRINT  "DISPLACEMENT  =  " ?NN1 

1430  LPRINT  "P  =  ";P 

1430  LPRINT  "w  =  "?A 

1530  LPRINT  "OMEGA  =  ";B 

1535  LPRINT  "VERT  SEPARATION  =  " ; Z 

1510  LPRINT  "NOISE  =  "?S1 

1515  LPRINT  "MAX  LATERAL  RANGE  =  "?RS 

1520  LPRINT  "NUMBER  OF  INCREMENTS  =  ";N7 

1525  LPRINT  ".LPRINT  "LTR  RNG    PD<CC)      PD<:SL)      PDCOPTV 

1533  L3=-R8 

1535  FOR  1=0  TO  NS 

1540  LPRINT  L2?"     "!Dl(I)i"      "?D2CI);M      U?D5<I) 

1545  L2=L2*D4 

1550  NEXT  I 

1560  GOTO  1050 

1605  PRINT  "END" 

1610  END 

1620  U=X0:V=Z 

1533  GOSUB  1200 

1640  D=J:H0=K 

1542  RH=«0. 1*P/<ORF*S2))*0. 333 

1543  SIG=AL*«H 

1545  X=< RH-H0)/SIG 

1546  GOSUB  1390 
1647  D5CE)=Y 

165a  3a=C0S<:B)*C0S(C3-A)  :J0=COS<:D)*CaS<:B)*SIN<C0-A>-SIN<:D)*SIN<B) 

1663  Na»-<siN<:D)>*cascB)*siN<;ca-A)  )-cos<:d:>*sin(B> 

1573  Bi«cas<F)*cas<:c0)  :ji=cos<d)*cos<f^*sin«:c3)-sin<:d)*sin«:f-) 

16SB  Nl=-<SIN<:D)*CaS(F)*SIN<:C3)  )-COS<.D)*SIN<F) 

1630  K'1=P/10/H0a3 

1700   A2=2*80*Bl-J3*Jl-N0*Nl:Al=3:*<N0>t-Bl-t-B0*Nl)  : A0=2*N0*N1-90*B1-J0*J1 

1710   S0=0 

1723   FOR    1=0    TO    M-l 

1730   S=«I-<M-l)/2)*03 

1740   Q=S/H0 
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1750  G=l/<:i+Q*GD  *2.  5 

1760  G=<A2*Q*Q+A1*Q+A0:>*G 

1770  G(I)=G:S0=S0+G*G 

1730  NEXT  I 

1730  RETURN 

1300  K2=SQR(S0) 

1310  K(:E.1=K2 

1320  VS=-V6+K1*SGR<S0)/S1 

1330  L0=K1*K1*S0/<S1*S1)  :G3=M+L0:B3=l-H_0/(ri+L0) 

1840  V3=-SQR<2*V7/B3)+SQR(2*A3/B3-1) 

1350  X=V3:G0SUB  1930 

1SE0  D1CE)»Y 

1370  X=V9:G0SUB  1330 

1330  D2<:E)=Y 

1330  RETURN 

1300  K=SQRiU*U+V*V) :IF  K=0  THEN  J=0:RETURN 

1905  UK=U/K:VK=V/K 

1907  IF  UK> 0.999999  AND  VK> 0. 399999  THEN  J=0: RETURN 

1903  IF  UK>0. 999999  THEN  J=-ATNCVK/SQR  <>VK*VK+1 ) )+90: RETURN 

1909  IF  VK> 0.999999  THEN  J=0: RETURN 

1910  MM=ATN<:UK/3QR<:-UK*UK+D)  :J=-ATN(VK/3QR<:-VK*VK+l))+90:  IF  MM<0  THEN  J=3E0-J 
1915  RETURN 

1920  Y=X:IF  X)0. 5  THEN  Y=l-Y 

1930  Y=3QR<L0GC1/Y/Y)) 

1940  G0-2. 515517:G1=0. 302S53:G2=0. 010323 

1950  HI =1. 432733 :H2=0. 139269 :H3=1. 30SE-03 

1960  Y=Y-<G0+Y*<  G1+G2*Y> > / < 1+Y*<H1+Y* < H2+H3+Y) ) ) 

1970  IF  X>0.5  THEN  Y=-Y 

I960  RETURN 

1990  Y=X:IF  X <0  THEN  Y=-Y 

2000  UI=l/<:i+0.231G419*Y) 

2010  Q1=0.  3133S153:Q2=-0. 356563732:03=1. 731 47793 :Q4=-1. 32 122553 :Q5=1. 33027442 

2020  IF  Y>24. 23  THEN  Y=0:GGTG  2070 

2025  PI=3. 14159265 

2030  Y=EXPC-<Y*Y/2)  ) /SQR<2*PI  )*W*CQ1+W*<  Q2+W*<:Q3+W*<Q4+W*Q5)  ) )  ) 

2070  IF  X>0  THEN  Y=l-Y 

2075  Y=C INT C 10000*Y) )/ 10000 

2030  RETURN 
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MAX  FREQ  -  1.3 

ADJ  INT  TIME  =  20.55555535 

SAMPLE  SIZE  =  37 

F/A  RATE  =  3 

PF  =  0.316S665SS6 

LAT  =  30 

LON  =  60 

PHI  =  59.43075979 

DIPOLE  COURSE  =  45 

DIPQLE  SPEED  =  10 

SENSOR  COURSE  =  315 

SENSOR  SPEED  =  220 

REL  COURSE  =  312.397439 

REL  SPEED  =220.227153 

EARTH  FIELD  =  52506.5513 

LONG  MOMENT  =  0 

TRAN  MOMENT  =  0 

VERT  MOMENT  =  0 

DISPLACEMENT  =  4000 

P  =  634634292 

w  =  32.53750751 

OMEGA  =  27. 11173179 

VERT  SEPARATION  =  200 

NOISE  =  3. 1 

MAX  LATERAL  RANGE  =  1500 

NUMBER  OF  INCREMENTS  =  15 


LTR  RNG 

PDCCC) 

PD<:SL)      PDC.OPT) 

-1500 

0. 0544 

0.0152      0 

-1400 

0. 0696 

0.0166      0 

-1300 

0. 0942 

0.0131      0 

-1200 

0. 1391 

0. 0209      0 

-1100 

0.2212 

0. 3255      3 

-1000 

0. 3755 

0. 0403      0 

-900 

0. 6399 

3.3831      0 

-S00 

0. 9191 

0. 2233      0 

-700 

0.999 

0. 5743      0 

-600 

1      0. 

9973      2E-34 

-500 

3.3133 

-400 

0. 1605 

-300 

0.567 

-200 

0.  2252 

-100 

0.9774 

0     1 

1 

0. 9292 

100     i 

L      1 

0. 9774 

200     1 

L      1 

0.2262 

300     3 

L      1 

0.567 

400    : 

L      1 

3. 1505 

500     3 

L      0.9999      0.0133 

600     G 

3.  9997 

0. 7654      2E-34 

700     C 

3.9157 

0.  2129      0 

200     ( 

3.  5775 

0. 0677      3 

900     ( 

3.3014 

3.  0331      3 

1000 

0. 1644 

3. 3225      0 

1100 

0. 1003 

0.3124      3 

1200 

0. 0523 

3.3165      3 

1300 

0. 050S 

3.3156      3 

1400 

0. 0405 

3.3151      3 

1500 

3.334 

0.0142      0 
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